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Three-dimensional iV-vector spin models may define universality 
classes for such diverse phenomena as i) the superfluid transition in 
liquid helium (currently investigated in the micro-gravity environment 
of the Space Shuttle) and ii) the transition from hadronic matter to a 
quark-gluon plasma, studied in heavy-ion collisions at the laboratories 
of Brookhaven and CERN. The models have been extensively studied 
both by field-theoretical and by statistical mechanical methods, includ- 
ing Monte Carlo simulations using cluster algorithms. These algorithms 
are applicable also in the presence of a magnetic field. Key quantities 
for the description of the transitions above — such as universal critical 
amplitude ratios and the location of the so-called pseudo-critical line 
— can be obtained from the models' magnetic equation of state, which 
relates magnetization, external magnetic field and temperature. Here 
we present an improved parametrization for the equation of state of 
the models, allowing a better fit to the numerical data. Our proposed 
form is inspired by perturbation theory, with coefficients determined 
nonperturbatively from fits to the data. 



I. INTRODUCTION 



The TV-vector (continuous-spin) models are proposed as class representatives for 
phase transitions in several interesting physical systems, such as the superfluid 
transition in liquid helium, in the N = 2 case [1], and the deconfinement transition 
in quantum chromodynamics (QCD) with 2 flavors of light quarks, in the N = 4 
case. The deconfinement phase transition is obtained when hadrons (e.g. protons 
and neutrons) melt into a quark-gluon plasma at very high temperatures, such as 
the temperatures that were present at beginning of the universe. There is great 
interest in describing this transition and in obtaining the properties of the high- 
temperature phase, a new state of matter that might be present today in the interior 
of neutron stars. In the case of two degenerate light quark flavors (i.e. up and 
down), the transition is believed to be described by the three-dimensional 4- vector 
model. More precisely, one invokes the effective a- model [2], a (three-dimensional) 
Ginzburg-Landau effective theory written assuming universality and respecting the 
chiral symmetry of QCD. The theory relates the QCD order parameter, which is 
the chiral condensate < tpxjj > (where tp is the quark field), to the magnetization 
of a continuous-spin model. The analogue of the magnetic field H is given by the 
(nonzero but small) quark mass and the reduced temperature is defined for lattice 
QCD as t ~ 6/<? 2 — 6/<7c 2 (0), where g is the lattice coupling. For two quark flavors 
one then obtains — if transition is second order — a three-dimensional 4-vector 
model in the presence of a magnetic field. The equivalence just described allows 
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one to study critical properties of the QCD phase transition from the spin-model 
equation of state. One of these properties is the so-called pseudo-critical line, the 
analogue of the critical point for the case of nonzero magnetic field. The pseudo- 
critical line is defined by the points where the susceptibility shows a (finite) peak, 
corresponding to the rounding of the divergence observed for H = 0, T = T c . 

We note that the equivalence between 2-flavor lattice QCD and the 4- vector model 
is still not verified in comparisons of the respective numerical data (see e.g. [3,4] 
and references therein). Thus, a better knowledge of the magnetic equation of state 
for the 4-vector model is of great importance to achieve higher precision in this 
comparison, to verify if the equivalence really holds and/or to establish the nature 
of the QCD phase transition, recently claimed to be of first order [5]. 

Also in the case of the 2- vector (or XY) model, a high-precision nonperturbative 
determination of the equation of state is of interest, since there are still discrepancies 
between the latest experimental and perturbative-renormalization-group values for 
critical quantities at the superfiuid helium transition [6]. Note that both these 
determinations are very accurate, while the available nonperturbative values (from 
Monte Carlo simulations at zero magnetic field [7,8]) are not as precise. 

The Hamiltonian for the TV-vector models is given by 

n = - J E s *- s s - H E s - w 

<i,j> i 

where the spin variables Si are taken as vectors on a sphere of unit radius in an 
./V-dimensional space. The main difference with respect to the Ising case is the pos- 
sibility of configurations where the spins are locally aligned but for long distances 
this alignment is lost, yielding a null average for the magnetization. Such configu- 
rations — called spin waves — possess arbitrarily low energy and tend to destroy 
the order of the system even at low temperatures. In d = 3 the models display a 
phase transition, with the presence of spontaneous magnetization below the criti- 
cal temperature, but the spin waves lead to Goldstone-mode induced singularities, 
causing a diverging susceptibility and strong finite-size effects for all T < T c when 
H->0. 

Monte Carlo simulations can be performed very efficiently for ./V-vector models, 
via cluster algorithms, also in a magnetic field [9-11]. The simulation in the presence 
of external field has the advantage that one can obtain the equation of state directly 
from the data, as described below. Also, in this case one can "measure" the actual 
magnetization of the system, without the need of an estimator such as the absolute 
value. In fact, the magnetic field already selects only one of the equivalent (zero- 
field) ground states that would lead to the same value of the estimator but would 
average to zero over the simulation in the zero-field case. The numerical simulation 
is done via cluster algorithm, which can be applied to the case of nonzero magnetic 
field by employing the ghost-spin trick. The observables are the magnetization 
parallel to the magnetic field and the two susceptibilities (parallel and orthogonal 
to H). 

We have recently proposed [12] an improved parametric form for the equation of 
state of the models. Our proposed form — inspired by perturbation theory — is a 
series expansion with two sets of terms, which contribute (mainly) separately to the 
description of the high- and low-temperature regions of the phase diagram. In this 
way we achieve a better description of the low-temperature phase at zero magnetic 
field (i.e. the coexistence line), characterized by the singularities described above. As 
a consequence, we are able to obtain a very precise characterization of the pseudo- 
critical line for the 4-vector model. We are currently applying the parametrization 
to a study of the N = 2 case. This will allow a better determination of the ratio of 
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critical amplitudes for the specific heat in the superfluid helium transition. Here we 
present preliminary results of this study, comparing the equations of state obtained 
for the cases N = 2, 4. 

II. THE MAGNETIC EQUATION OF STATE 

At infinite volume, the scaling function for the singular part of the free energy is 
given by 

F s (t,h) = b~ d F s {b y H,b yh h) , (2) 

where b is arbitrary, t = (T — T c )/Tq, h = H/Hq and yt, yh are related to the 
usual critical exponents j3 and 7. The above form implies the relation between 
magnetization and the applied magnetic field, known as the magnetic equation of 
state 

M/h 1 '* = f M {t/hW) . (3) 

Equivalently 

y = fix) , (4) 

with y = h/M s , x = t 13 /M. The normalization constants are given by /(0) = 1, 
/(— 1) = 0. In the case of the N- vector models, the singularities at low temperature 
determine the behavior of the magnetization as the square root of H. This behavior 
(i.e. the Goldstone-mode effect) is included in the following Ansatz [13] for the 
equation of state at low values of x 

x = -l + ay 1/2 + by + cy 3/2 + ■■■ (5) 

Note that these effects are present in -/V-vector models along the coexistence line, 
i.e. at low temperatures and small magnetic field (or equivalently, at low values of 
the variable a;). 

The pseudo-critical line, described in the previous section, is given by finite peaks 
in the susceptibility. It characterizes the critical region when the external field is 
not zero (e.g. in the QCD case). The scaling form for the susceptibility along the 
pseudo-critical line is given by 

X = dM/dH = (1/Ho) h 1 ' 5 - 1 f x (t/hW) . (6) 

Note the \ has a peak at t p for each fixed h and that 

tp^zph 1 ^ 5 (7) 
H 0Xp = h 1 ' s - 1 f x (z p ). (8) 

Thus, the location of pseudo-critical line is given by z p , obtained from the scaling 
function (or equation of state) for the susceptibility, which involves the derivative 
of fiviit/h 1 ^ 5 ). Note also that z p is a universal constant. 

For the 4- vector and the 2-vector models, the equation of state was determined 
numerically respectively in [9] and [10,11], by taking into account the Goldstone- 
mode singularities and determining the location of the pseudo-critical line. The 
fitting function for f(x) = y used an interpolation of two forms: the Ansatz in Eq. 
(5) at low x and Griffiths's analyticity conditions at large x 

x{y) = AyV~< + By^-WIt +■■■ (9) 
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Then the equation of state fu{z) = M/h}/ 5 is obtained from x(y). The problem 
is that the "transition" between the two fitting forms above is close to the pseudo- 
critical point, which is itself obtained from the equation of state for the susceptibility 
[defined in Eq. (6) above], involving a derivative. It would be therefore preferable 
to use a parametric form for f(x) without the need for interpolation, ensuring a 
smooth derivative and a better determination of the pseudo-critical point z p . To 
this end we can consider the parametric form 

M = m {) R t) 6, (10) 

t = r (i -e 2 ) , (ii) 

H = h R fjS h{9) . (12) 

This yields 

where h{9) is an odd function, with root given by 9q. This form was introduced in 
[14] by Guida & Zinn- Justin for the Ising model (therefore without considering the 
effect due to Goldstone modes at low x) and used for perturbative studies of the 
TV-vector case with the inclusion of the leading contribution of 9q. The form leads 
to a smooth curve, and allows a direct relationship with critical amplitude ratios. 
We propose [12] an improved parametrization given by 

n m 

h(9) = 9(l- 9 2 /9 2 ) 2 (l + J2^) x[l + J2d 3 (l- 9 2 /9 2 y}. (14) 

i=l j=l 

The above form is based on the parametrization used perturbatively in [14] for the 
Ising model, but takes into account terms associated with the effects of singularities 
induced by Goldstone modes, as discussed in the Introduction. These new terms 
are included by means of the dj coefficients, associated with an expansion around 
the coexistence line. (The dj's are considered in addition to the usual Cj coefficients, 
related to the high-temperature/high- x behavior.) 

In Ref. [12] we have used the proposed form above for fits to existing Monte Carlo 
data for the 4- vector case. From our fits we see that dj's are indeed necessary for 
the description of the data. Our best fit is obtained considering (in addition to 6>o) 
two coefficients of type c and two of type d 

9 2 a = 2.17(4) (15) 
ci = 0.9(1), c 2 = -0.62(7) (16) 
di = -1.56(4), d 2 = 1.15(5). (17) 

We note that these results have errors that are one order of magnitude smaller than 
the perturbative description, and fit to the data with a value of \ 2 P er degree of 
freedom that is two orders of magnitude smaller. The fit is shown together with the 
data in Fig. 1 (left panel). 

As explained above, the location of the pseudo-critical line (useful for comparison 
between QCD data and the iV-vector model's equation of state) is obtained from 
the peak in the scaling function for the susceptibility (see Eq. 6). The peak can 
be determined numerically from the two equations (13) by varying 9. Our result is 
obtained with less than 1% of error 

9 P = 0.587(2), z p = 1.29(1), f x (z p ) = 0.341(1). (18) 
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The results are in agreement with previous determinations of z p and f x (z p ), but 
our error for z p is much smaller. 

Next, we show our preliminary results for the N — 2 case. We have produced 
new data, simulating the three-dimensional ATF-model in the presence of a magnetic 
field by means of the (Wolff) cluster algorithm. We use the improved form of the 
model's Hamiltonian, introduced by Hasenbusch and Torok for the zero-field case 
in [7]. We also use their values for the critical temperature and critical exponents. 
When fitting the data to our parametrization of the equation of state, we find 
— regarding the role played by the c and d coefficients — essentially the same 
characteristics as in the N — 4 case described above, with the difference that in this 
case 6 parameters are needed. Our best fit is obtained using 

6l = 3.25(2) (19) 
ci = 1.05(4), c 2 = -0.11(3), c 3 = 0.53(2) (20) 
di = -6.75(2), d 2 = 14.7(2). (21) 

We show the data together with the fit in Fig. 1 (right panel). We see that the slope 
of the curve is significantly higher for the N = 4 case, corresponding to a stronger 
effect of Goldstone-mode singularities, as has already been found in [10,11]. Notice 
that the data in the N — 2 case have smaller error bars, leading to a very precise 
determination of the curve. This will enable us to calculate (see e.g. [12]) the critical 
amplitude ratio of the specific heat with the same accuracy as the experimental and 
perturbative values. We quote our values for the location of the pseudo-critical 
region in the 2- vector case 

9 P = 0.563(2), z p = 1.61(2), f x (z p ) = 0.349(1). (22) 

As for the A" = 4 case, the results are in agreement with previous determinations 
of the pseudo-critical line [10,11], but the error for z v is smaller by one order of 
magnitude. 




FIG. 1. Plot of the data together with the fitting form for y(x) in the 4-vector case (left 
panel) and in the 2-vector case (right panel). 



III. CONCLUSIONS 

We have introduced an improved parametric form for the description of the equa- 
tion of state of 3d A'-vector models. We show that the new parametric form indeed 



5 



provides a better fit to the numerical data as compared to previous parametriza- 
tions. In particular, the consideration of the dj coefficients is essential for a good 
description of the Monte Carlo data in the whole range of values of x. Also, we 
were able to verify clearly the different roles played by a and dj parameters in the 
high- and low- a; regions. Indeed, in this form the coefficients c, and dj contribute 
respectively to the high- (9 « 0) and low- (6 w 6q) temperature regions. We also 
stress that, in addition to providing a better fit to the numerical data, the expres- 
sion considered is a continuous function, needing no interpolation between the two 
x regions. This is particularly useful for the determination of the pseudo-critical 
line, since the interpolating form introduced in [9] is unstable precisely in this re- 
gion. In fact, our determination of z v is very precise in comparison to the previous 
estimates from the interpolated form and the perturbative equation of state. As a 
consequence of a better determination of the pseudo-critical line in the N = 4 case, 
one may get an unambiguous normalization of QCD data for comparison to the 
4- vector equation of state, showing better agreement for larger quark masses [3,4]. 
We are currently extending our analysis to the N = 2 case. Our preliminary results, 
presented here, confirm the fact that the equation of state can be obtained with very 
high precision using our method. This will allow a determination of the the critical 
amplitude ratio of the specific heat with the same accuracy as the experimental and 
perturbative values, which is of a few tenths of a percent [6]. 
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